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Abstract: This chapter defines and implements a non-Bayesian 
fusion rule for combining densities of probabilities, derived from im- 
precise knowledge. This rule is the restriction to a strict probabilis- 
tic paradigm of the Proportional Conflict Redistribution rule no 6 
(PCR6) developed in the DSmT framework for fusing basic belief as- 
signments. A sampling method for probabilistic PCR6 (p-PCR6) is 
defined. It is shown that p-PCR6 allows to keep the modes of local 
densities and preserve as much as possible the whole information in- 
herent to each densities to combine. In particular, p-PCR6 is able 
of maintaining multiple hypotheses/modes, when they are too distant 
for fusion, contrariwise to classical technique. The question of se- 
quential filtering by p-PC'R6 is addressed, thus implying the necessity 
to handle the redundancy of the information. 
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Notations 


e |x = y| is the Dirac distribution of variable x on value y, 


e [|b|, function of Boolean b, is defined by [true] = 1 and /{false| = 0. In 
particular, J|r = y] could be seen as a discrete counterpart of the Dirac 
opes. 


4.1 Introduction 


Bayesian inference is a powerful principle for modeling and manipulating prob- 
abilistic information. In many cases, Bayesian inference is considered as an 
optimal and legitimate rule for inferring such information. Bayesian filters for 
example, and their approximations by means of sequential Monte-Carlo, are 
typically regarded as optimal filters |1, 2, 7]. 


However, Bayesian methods need strong hypotheses, in particular about 
the information prior and the independence prior. A degradation of the per- 
formance of Bayesian filter occurs if the filter is not correctly initialized or 
updated, in accordance to the models in use. Being given a model of the sys- 
tem kinematic and of the measurement process, the main issue is to develop 
filtering methods which are sufficiently robust against the bias at the initializa- 
tion as well as error in modeling. In this paper, a non-Bayesian rule for fusing 
the probabilistic information is proposed. This rule, denoted p-PCR6, is the 
restriction to the probabilistic paradigm of the Proportional Conflict redistri- 
bution rule no.6 (PCR6) which has been proposed in [12] for combining basic 
belief assignments. p-PCR6 is also an extension of discrete PCR6 version to 
its continuous probabilistic counterpart. 








PCR6 has been first established for combining evidences (i.e. discrete be- 
lief assignments) in the DSmT framework. In particular, it has been designed 
in order to cope with highly conflicting and uncertain information. This rule 
could be considered in a probabilistic paradigm by restricting the basic be- 
lief assignment involved to only probabilistic belief assignment! , and directly 
extended to densities of probabilities. This rule is non-Bayesian by nature. Al- 
though Bayesian techniques are widely well known and used in target tracking 
community (including authors works in tracking), it is interesting to see how 
such new approach can perform to estimate its real interest and potentiality. 


'The denomination probabilistic belief assignment is preferred to Bayesian belief as- 
signment, generally used in the literature, since we consider that Probability and Bayesian 
inference are distinguishable notions. 
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Surprisingly, it turns out through our works, that such approach is robust to an 
erroneous modeling: in particular, it is able of maintaining multiple hypothe- 
ses, when they are too distant for fusion. The resulting p-PCR6-based filter 
happens to be essentially non-linear, and has been implemented in our simula- 
tion using particle filtering techniques. In particular, the p-PCR6 multisensor 
filter developed here is based on a quite simple and direct implementation in 
terms of particles drawing and resampling. At the end of this chapter, the 
question of sequential filtering is addressed. In this case, it is necessary to take 
into account the redundancy of information over the time. Then, p-PCR6 is 
adapted in order to remove this redundancy. 





Section 4.2 introduces the PCR6 rules, and establishes some results about 
probabilistic PCR6. A sampling method is deduced. Section 4.3 compares the 
results of the Bayesian rule and of probabilistic PCR6 on a simple example. On 
the basis of this comparison, some arguments about the robustness of PCR6 
are given. Section 4.4 investigates the sequential filtering issue. Application of 
p-PCR6 to distributed filtering is provided as example. Section 4.5 concludes. 


4.2 PCR6 formula for probabilities 


4.2.1 Definition and justification of PCR6 


The Proportional Conflict Redistribution rule no. 6 (PCR6) of combination [5] 
is an extension of rule PCR5 |10, 11]. These rules come from the necessity 
to manage precisely and efficiently the partial conflicts when combining con- 
flicting and uncertain information expressed in terms of (quantitative) belief 
assignments. These rule have been proved useful and powerful in several ap- 
plications where it has been used [12]. PCR5 and PCR6 are equivalent, when 
restricted to only two sources of information. 





Let be given an universe of events O. A distribution of evidence over O 
is characterized by means of a basic belief assignment (bba) m : P(O) — R+ 
such that: 


m(@)=0 and » Ix eI. 


Xco 


where P(O) is the set of subset of O.? 


2In the general case, bba could also be defined over hyper-power sets (Dedekind's lat- 
tice) [12]. 
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A bba typically represents the knowledge, which can be both uncertain and 
imprecise, that a sensor provides about its belief in the true state of the uni- 
verse. The question then arising is How to fuse the bba’s related to multiple 
sensor responses? The main idea is to corroborate the information of each 
sensor in a conjunctive way. 


Example: Let's assume two sources with basic belief assignments m1 and mo 
such that m4(A) = 0.6, mı (AU B) = 0.4 and m2(B) = 0.3, m2(AU B) = 0.7. 
The fused bba is then characterized in a conjunctive way by: 


ma (AN B) = mi(A)mo(B) = 0.18, 
ma(A 2) 7 mi(A Ayno(AU B) = 0.42, 
ue m1 (AU B}mo(AU B) = 028 . 


The conjunctive consensus works well when there is no possibility of con- 
flict. Now, make the hypothesis AN B = Ø. Then, it is obtained mA(()) = 0.18, 
which is not an acceptable result for a conventional interpretation of Ü as a 
contradiction. Most existing rules solve this issue by redistributing the con- 
flict ma (Ø) over the other propositions. In PCR6, the partial conflicting mass 
m1(A)m»(B) is redistributed to A and B only with the respective proportions 
xa = 0.12 and rp = 0.06, according to the proportionalization principle: 

XA LB m1(A)m2(B) 0.18 


mi(A) E m»(B) i m1(AÀ) + m2(B) 7 09 — "e 





Basically, the idea of PCR6 is to transfer the conflicting mass only to the ele- 
ments involved in the conflict and proportionally to their individual masses. 


Some theoretical considerations and justifications already briefly aforemen- 
tioned led to the following PCR6 combination rule. Being given two bba's 
m and mo», the fused bba mpong according to PCR6, or equivalently to PCR5 
in this case, is defined for any X € P(O) \ {Ø} by: 


MPCR5/PCR6 (X) = MA (X) 


where m,(-) corresponds to the conjunctive consensus: 


ma(X) 2 m1(Y1)m2(Y2) . 
YNY=X 
Y1,Y2EP(O) 
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When fusing s > 2 sources of informations, characterized by the bba’s m1 to 
Ms, the fused bba is defined in [5] by: 


MPCR6 (X) = MA (X) 


s—1 
]H meo (o) 
j=l 


+ 2s mi(X)? > p — , (4.2) 


DkciYego0X-9 | m;(X) + mot 
ioi. ePi i ) 2. T G0) 


where m,(-) corresponds to the conjunctive consensus: 


m(X)= M  [[m(Y2, 


Y1N -QY =X 21 
Xp. ,Y:€P(6) 


and the function o; counts from 1 to s avoiding 1: 


oj) =F x MF <Y+ Gt) xed. 


4.2.2 Reformulation of PCR6 


Definition (4.2) could be reformulated into a more intuitive expression. 


i=1 s _Y.=0 (V. 
"s cP(©) 2- mj (Yj) 
and then: 


mpore(X)=ma(X)+ SY. [[miy) =. (43) 


s Y =Ø :—1 | | 
en 2 mj) 
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At last, a new formulation of PCR6 is derived for X € P(O) \ (0): 


mpcne(X) = >, Tes: ») Fpecne(X ]|Yi: D , 


Yi, ,Y:€P(6) 


where the function Fpcre is defined by: 


Fpors (X |Y1:s) = 1 


Xt 


a ———— —— 7 





ny - 























When considering probabilistic densities instead of belief functions, the com- 
ponents [];_, m;(Y;) and jani RIDE RON) mi) have a 
straightforward interpretation. The first is interpreted as an independent gen- 
eration of answers by each source of information. The second is interpreted as a 
random choice among the answers or the consensus, weighted by the respective 
evidences. 


4.2.3 Definition of probabilistic PCR6 (p-PCR6) 


In [12], Dezert and Smarandache proposed a probabilistic version of the PCR5 / 
PCR6 rule (4.1) for two sources, by restricting the bba’s m1 and m2 to discrete 
probabilities P, and P5 which are called then probabilistic belief assignments 
(or masse’). Probabilistic belief masses are bba’s, which focal elements? con- 
sist only in elements of the frame O, i.e. the singletons only. When dealing 
with probabilistic belief assignments mı = P, and mə = P5, the conjunctive 
consensus is restricted to the same singleton, so that mA(X) = P,(X)P2(X). 


?Focal elements are elements of P(O) having a strictly positive mass. 
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As a consequence, the PCR5/PCR6 formula (4.1) for two sources reduces to: 
PA(X)PX(Y) 
P, X) = Pi(X)Po(X) + P(X —— 
ponsen (X) = PORO KO D grep 
YEO\{xX } 
P(X) P(Y 
+P(X M aA d | 
YeEo\{x} 
Now, it happens that: 
Pi(X) P(X) 
Ep PEDE 


P2(X)P,(X) 


P(X) P(X) = P(X) Pi(X) + Px(X) ’ 


+ P(X) 
and finally: 


Pecrs/pcre(X X) x = 72 9 7 ne ROSA 
zzu i 
(4.5) 


Of course, this formula generalizes in the case of PCR6 for any number of 
sources. Since: 


and, Tor X, Y, qj € 9, 


s—1 
( | You) NX Z0 if and only if Y; = — Yo =X, (4.6) 
k=1 
it comes: 
S lI E Vois) 
j—1 
Prons(X) = 9 P(X} >, [0 nn 
i=1 Far sYote=4) SO P; (X) + >. Lo aa 
j=l 


(4.7) 


Equations (4.5) and (4.7) are however difficult to handle practically. The re- 
formulated definition of p-PCR6 is introduced now. 


144 Chapter 4: Probabilistic PCR6 fusion rule 


4.2.4 Reformulation of p-PCR6 
From (4.4), it is deduced: 





Pecre(X) — 
. (uon = = Y,] +I Av =o x =y) DO 
= k=1 
D [Ro = 
Yi,,Y,cO i-l 2 PY) 
j=l 
Now, the property (4.6) implies: 
IX =Y¥, =---=Y.) +1] (| Ye =0| x =Yj] =1[X 2 Yj]. 
k=1 








As a consequence, the p-PCR6 rules is equivalently defined by: 


Ppcre (X) E > [AY = 
Yi, ,Y;cO i=1 » PO) 


(4.8) 


4.2.5 Extension of p-PCR6 on continuous propositions 


The previous discrete p-PCR6 formula is now extended to densities of proba- 
bilities of random variables. Formula (4.7) is thus adapted for the fusion of 
continuous densities p4,:-- , ps: 


Se 
pile) |] Poi Won) ,., 
= 
prosa) E ^p) | ——ÀÀÀ— — TJ due (49 
pe 


s—1 


LI 
pi(z) + 9 pono (Yor)? 
j=l 


Notice that p;(x) is put inside the integration, so as to deal with possible 
singularities, when p;(x) = 0. It is also necessary to prove that ppcne is a 
probabilistic density. And of course, it is possible to guess a reformulated 
definition of p-PCR6 for densities by means of (4.8). But, we establish now 
these results by calculus. First at all, a result is proved for computing the 
expectation based on ppocnse. 
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4.2.5.1 Expectation 


The expectation of a function according to the fused probability ppong is ex- 
pressed from the initial probabilities p1,---: , Ps by: 


S 


s—1 
pily) lH »so (9o) s—1 


=> | oto f — = fty,z) lH dvo dy 


jex] 
pily) + 9 Pog Uog) i 
j=l 


OOO 


Corollary. The density ppcre is actually probabilistic, since it is derived: 


| Prcré() dy = 1, 
o 


by taking f — 1. 
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4.2.5.2 Reformulated definition 


Prore(Z) = L (Tiro) npn.) | [ du, 


» pi(yi)ôly: =g] (4.11) 


where r(2ly15) = = 


Proof. 
Apply lemma 1 to the Dirac distribution f(y, z) = ó|y = z]. 
OOL 


4.2.6 Sampling method 


Being able to sample p;,--- , Ps, it is possible to sample ppcre by applying 
the definition (4.11). The implied sampling process (let z be the sample to be 
generated) is sketched as follows: 


1. For any k € {1,---,s}, generate y; according to pz, together with its 
evaluation pi (y), 


2. Generate 0 € [0,1] according to the uniform law, 


X PRO PRY k) = Xi P) 1 Dx (Yk) 


« 0 


3. Find 7 such that 
pat Pk(Yk) os | Pk(Yk) | 


4. Set z = yj. 





It is seen subsequently that p-PCR6 does not preserve the Gaussian distribu- 
tions. As a consequence, its manipulation is essentially addressed by means of 
a Monte-Carlo method, and the previous sampling method is implemented in 
the applications. 


The next section is devoted to a comparison of p-PCR6 and Bayesian rules on 
very simple examples. 
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4.3 Bayes versus p-PCR6 


4.3.1 Bayesian fusion rule 





In this section, we are interested in the fusion of two independent estimators 
by means of the Bayesian inference. Such fusion has to take into account the 
prior about the state of the system. Subsequently, this prior is chosen to be 
uniform. Although this is just a particular case of application, it is sufficient for 
our purpose, 2.e. the illustration of essential differences between the Bayesian 
and PCR6 approaches. 





4.3.1.1 General case 


In Bayesian filter, the estimator is explained by means of the posterior pro- 
bability p(x|z1,22) conditionally to the observation 21 and z2. Notice that 
this posterior estimation should not be confounded with the true state of the 
system. Now, our purpose here is to derive a rule for deriving the global es- 
timator p(x|21,22) from the partial estimators p(x|z1) and p(a|z2). Applying 
Bayes’ rule, one gets p(x|21,22) x p(z1, 22|r)p(x) .* To go further in the deriva- 
tion, one must assume the conditional independence between the two proba- 
bilistic sources/densities, i.e. p(z1,zo|x) = p(z1|v)p(z2|x). As a consequence, 
p(z|z1, 22) e p(z1|z)p(22|v)p(x) , and then: 


p(x|z1)p(x|22) 

p(z) | 
So, in order to compute p(x|21,22), it is needed both p(x|z1), p(a|z2) and the 
prior p(x). If one assumes uniform prior for p(x), and using notations Ppayes = 
p(-|21, 22), pı = plz) and po = p(-|z2), the Bayes’ fusion formula (4.12) 
becomes: 





p(x|z1, z2) x (4.12) 


De (0) OX pi(x)po(x) . (4.13) 
(It is noticed that a discrete counterpart of this result could also be obtained by 
applying Dempster Shafer rule to probabilistic belief masses) 


4.3.1.2 Gaussian subcase 





We investigate here the solution of the problem when pı and pə are Gaus- 
sian distributions. So assume for simplicity that p1(x) and po(r) are mono- 
dimensional Gaussian distributions: 


1 1 (x—31)? 4 1 1 (1—22)? 
T) = ——e 1 an p= 
n O1 V 2T Oo V 2T 


#p(a|B) x y means “p(a|B) is proportional to y for B fixed”. 
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In absence of prior information, one assumes p(x) uniform. The Bayesian rule 
requires to compute (4.13). Then, it is easily shown that Pgayes is Gaussian: 


= 2 
1 (T—TBayes) 
]  -iT Eee 


ayes \L) = — 
di > ) Obaya Van 


with 
2 oe 103 
Bayes c? ES c2 ? 
and E 
X — g? a E) 
Bayes Bayes (5 2) 


When oi = o2 = Ø, it is implied then: 


2 == = 
O B Ti F Tə 
eC — = and Tiye” ES 


Bayes 
Hence, the resulting standard deviation op, after Bayes fusion is equal to 
the initial standard deviation divided by the factor V2 and thus Opus €. 
This fusion process is optimal, when the model parameters are correct. Now, 
imagine that the difference zə — zı is obtained from a bias error of the model. 
For example, let us consider that the estimation of sensor 1 is correct but 
that the estimation of sensor 2 is erroneous, in regards to the deviation oc. 
Assuming x being the true state of the system, it comes most likely: pi(x) > 
PBayes(“@) > po(x). Thus, the Bayesian fusion propagates the errors. This 
implies an irrelevant estimation. It is noticed however, that the bias is divided 
by two, each time a fusion with a good estimation occurs, while the deviation 
is only divided by V2. Then, good estimations will make the process converge 
correctly after some iteration. 


The theoretical plots and those obtained with Monte Carlo simulation are 
given in figures 4.1, 4.2 and 4.3. These figures make the comparison with the p- 
PCR6 fused densities. This comparison will be discussed subsequently. It is yet 
confirmed that the Bayesian rule just concentrates the information, by reducing 
the deviation, even when the information are distant (that is putatively false). 
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Figure 4.1: p-PCR6 fusion versus Bayesian fusion (theoretical). 
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Figure 4.2: p-PCR6 fusion versus Bayesian fusion (theoretical). 
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Figure 4.3: p-PCR6 fusion versus Bayesian fusion (based on 10000 samples). 


4.3.2 Fusion based on p-PCR6 for Gaussian distributions 


The same Gaussian distribution, p, and po, are considered, but are now fused by 
p-PCR6 rule (4.9), thus resulting in density ppcno. The fused densities are both 
computed, figures 4.1 and 4.2, and sampled, figure 4.3. Direct computations 
are expensive, and are obtained in two steps: 





e Compute 1,(x) = Lo dy, where s € {1,2} and 5 € {1,2} V {s}, 


e Then compute ppecne(z) = pı (£) (x) + po(x)Io(x) . 


It appears clearly that computed and sampled densities match well, thus con- 
firming the rightness of our sampling method. Now, contrariwise to the Bayesian 
rule, it is noticed two different behaviors (which are foreseeable mathemati- 
cally): 


e When the densities p, and p» are close, ppone acts as an amplifier of the 
information by reducing the variance. However, this phenomena is weaker 
than for Ppayes. p-PCR6 is thus able to amplify the fused information, 
but is less powerful than the Bayesian rule in this task. 


e When the densities pı and p» are distant, ppcreg keeps both modes present 
in each density and preserves the richness of information by not merging 
both densities into only one (unimodal) Gaussian density. This is a very 
interesting and new property from a theoretical point of view, which 
presents advantages for practical applications. 
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In regards to these differences, it is thus foreseeable that the p-PCR6 should 
be more robust to potential errors. 


4.4 A distributed sequential filtering application 


4.4.1 Whitened p-PCR6 rule 


It has been seen that the p-PCR6 fusion of the same densities py = po will 
result in an amplified density ppcne. Of course, this is not practicable when 
the densities p; and p» are related to correlated variables. Consider for exam- 
ple that the state y are measured by z! and z?. The (distributed) posterior 
probabilities are ps(y) = p(y|z*) e p(y)p(z?|y) for s = 1,2. It happens that 
pı and p» are correlated, so that p-PCR6 should not be applied directly. In 
particular, the fusion of p; and po by means of p-PCR6 results in a density 
Dpore Stronger than the prior p over y, even when there is no informative mea- 
sure, i.e. p(z*|y) = p(z*)! In order to handle this difficulty, we propose a 
whitened p-PCR6 rule, producing a fused density Pwhitepcre from the updated 
information only: 


Pwhitepcre ly) = J L pi(yi)p2(y2)r(yly1, y2) dyidy , 


p(yi|z!) _ p(yo|z2) (4.14) 
? lary OT Gay Ôly2 = y 
where n(y|yi; y2) = plz") , plz?) i 
ply) +  ply2) | 


In (4.14), the proportion 2 Js should be considered as the information intrin- 





sically obtained from sensor s. It happens that the whitened p-PCR6 does not 
change the prior when there is no informative measure, i.e. Pwhitepcre(Y) = ply) 
when p(z?|y) = p(z*) for s = 1,2. 


4.4.2 Theoretical setting 


A target is moving according to a known Markov prior law. Let y; be the state 
of the target at time t. It is assumed: 


D(y1:t41) = Dia yop) . 


In order to estimate the state of the target, S sensors are providing some 
measurements. Denote z? the measurement of the state y; by sensor s. The 
measure is characterized by the law p(z?|y:), which is known. It is assumed 
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that the measure are made independently, conditionally to the state: 


S 


par lue) = [pa lu) - 


= | 


Our purpose is to derive or approximate the optimal estimator, p(yt+1 Ira 
from the distributed retroacted estimators, p(y1l21,z2f,1), related to sen- 
sors s. There is a Bayesian approach to this problem, and we propose some 
comparison with a p-PCR6 approach and a whitened p-PCR6 approach. 


4.4.2.1 Distributed Bayesian filter 
This filter is derived from: 


p(yisilzii)- J p(ysilu)p(ulzre ) due , (4.15) 
Yt 
1:9 s S 1:9 
p(yesailzii ru) CX p(zzaalyeci)p(yena zia ) D (4.16) 
( | 1:9 ) (I p(yetilZie, — i ( | d (4 i7) 
P\Ut+1|%1:t4+1) & : DVUt-4-1121: à ; 
B kun p(yexilzii) 


This approach is unstable, when some components of the target state are non- 
observable; for example, adaptations of the method are necessary |2| for bearing 
only sensors. However, the method will be applied as it is here to bearing only 
sensors, in order to compare to the robustness of the PCR6 approach. 


4.4.2.2 p-PCR6 filter 
This filter is derived from (4.15), (4.16) and: 


S 
(I P(Yi+1 ie si) A (Yeti luz) dy 
(4.18) 
Ya palet. 2241) 6lyet = vial 
3e P(y£ ale. ze41) 


1S 
P(ye+1(21:¢41) - | 
125. 
Yer 


where m(ye+ilyii) = 


1 


and p(y£,4|zE? , 22,4) is an instance of p(yz41|21:7, 22,1), obtained by just re- 
placing yr+1 by yey. 

It is noticed that this filter is necessary suboptimal, since it makes use of the p- 
PCR6 rule on correlated variables. The whitened p-PCR6 filter will resolve this 
difficulty. However, it is seen that the p- PCR filter still works experimentally 
on the considered examples. 
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4.4.2.5 Whitened p-PCR6 filter 
This filter is derived from (4.15), (4.16) and: 


S 
(II p (yide, a) re luli dii 


s—1 


1:S 
p(yt+1 21441) -[ 
Iss 
Yeti 


p» m PUE p3 la t41) 5 (4.19) 


— „S 
P? alzi?) cei = Yet] 


where m(ye4i|Yeti) = WET 
pe P(Uf ali pua) 
= palal) 


Again, yz,4, is just an instance of y;41 for sensor s . 


These filters have been implemented by means of particles. The sampling of p- 
PCR6 has been explained yet, but it is not the purpose of this paper to explain 
all the theory of particle filtering; a consultation of the literature, e.g. [9], is 
expected. 


4.4.3 Scenario and tests 





These examples are retrieved from |4]. This work has been implemented by 
Alois Kirchner during his internship in our team. 


4.4.3.1 Scenario for passive multi-sensor target tracking 


In order to test the p-PCR6 fusion rule, we simulate the following scenario: in 
a 2-dimensional space, two independent passive sensors are located in (0,100) 
and (100,0) in Cartesian coordinates. These sensors provide a noisy azimuth 
measurement (0.01 rad. normal noise) on the position of a moving target. We 
associate a tracking particle filter to each sensor. The motion model is the 
following : 

dud = Tt + 0.1 x N(0, 1) 

Yer. = ÿ+0.1xN(0,1) 

Tt+1 = £+ + dt x t+ + 0.3 x N (0,1) 

Yt+1 = Yt + dt gi + 0.3 x N(0, 1) 


where dt = 1 time unit and N(0,1) is the normal distribution. 


(4.20) 


In our simulations, each local particle filter is implemented by means of 200 
particles. At each time step, we proceed to the fusion of the local posterior densities 
and then re-inject the fused state density into each local filter (feedback loop). Three 
different paradigms are considered for the fusion: Bayesian, p-PCR6 and whitened 
p-PCR6 rules. These filters try to estimate both the position and speed of the target 
which is assumed to follow a quasi-constant velocity model. It is noticed that we 
are dealing directly with both the observable and non-observable components of the 
target state. 
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4.4.3.2 A simple example 





In this first example, the filters are well initialized (we give them good starting speed 
and position). The mobile follows a non-linear trajectory (figure 4.4), in order to show 
the capability of this distributed filter to converge. On this example, the Bayesian 
filter manages to track the target with some difficulties during the last curve in 
figure 4.4. On the same example, p-PCR6 and whitened p-PCR6 rules have been 
tested with success. While both filters have to reestimate the speed direction at each 
turn, it appears that this reestimation is more difficult for p-PCR6. This difference 
is also particularly apparent during the last curve. 








True trajectory 


Bayes 
Whitened p-PCRS ——— 
p-PCR5 —— — 





Figure 4.4: Estimated trajectories using different tracking methods. 


Figure 4.5 displays the particle cloud of the whitened PCR6 filter during and 
after the last curve. The variance rises along the curve, resulting in the cross-like 
cloud of sub-figure 4.5(a), which is characteristic to the p-PCR6 fusion: the branches 
correspond to the direction the sensors are looking at. Then, the p-PCR6, by ampli- 
tying the zone where the filters are according to see the object, allows the process to 
converge again toward the object real position in an expansion-contraction pattern 


(figure 4.5(b)). 
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True trajectory 
Whitened p-PCRS ——— 


Particles + 





(a) Time step 160. 


True trajectory 
Whitened p-PCRS ——— 


Particles + 





(b) Time step 170. 


Figure 4.5: Particle clouds for whitened p-PCR6 in the last curve. 
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In more difficult cases, with poor initialization for instance (see figure 4.6), both 
p-PCR6 and whitened p-PCR6 manage to follow the target, while the Bayesian filter 
diverges in about one third of the cases and give mitigate results otherwise. 


True trajectory 


Bayes 
Whitened p-PCRS ——— 
p-PCR5 —— 





Figure 4.6: Estimated trajectories using different tracking methods. Poor ini- 
tialization: null speed and 10 units away starting position. 


Next sections investigate more thoroughly the properties of the whitened p-PCR6 
filtering. 


4.4.3.3 Whitened p-PCR6 robustness against poor initialization 


In order to test the capability of (whitened) p-PCR6 to recover from erroneous mea- 
surements or from a total contradiction of the local estimations, we considered two 
scenarios in which the filters are badly initialized at various degrees. In these scenar- 
ios, the real trajectory of the object is the same: it starts from (200,0) and moves 
toward (200, 150) at a constant speed (0,1). 

In the first scenario (figure 4.7), the first filter, which sensor is placed at (0, 100), 
is initialized at position (190,10) and at speed (0,0). The second filter, which sensor 
is at (100, 0), is initialized at position (210, 10) and at the same speed (figure 4.7(a)). 
As the estimated positions are far from the real one and both sensors are looking 
at the object from a remote position, the particle cloud quickly spread horizontally 
(figure 4.7(b)). Then the (whitened) PCR6 begins to find zones where both filters es- 
timate a non-negligible probability of presence and amplifies them until convergence 
(figure 4.7(c)). Though the particle cloud still seems to be fairly spread (because of 
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sensors remote position), the global estimate is very close from the real position and 
speed, and will remain so until the last time step (figure 4.7(d)). 


— | x |y | x speed | y speed | 
First | Filter 1 | 190 | 10 | 0 | 0 — 


example [ Fier? 210 [10] o | o 
example 


Table 4.1: Initialization data. 





True trajectory —— 
Whitened p-PCR5 — — 
Particles + 





(a) Time step 1. (b) Time step 10. 








(c) Time step 20. (d) Time step 60. 


Figure 4.7: The real mobile starts at (200,0) and moves upward at constant 
speed (0, 1); poor filters initialization. 
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Our second example (figure 4.8) is an extreme case: the initialization is quite worse 
(see table 4.1), since our motion model assumes nearly constant speed and therefore 
makes it hard to recover from such erroneous and contradictory speed initialization. 
An interesting point is that, for a tight prediction noise, p-PCR6 sometimes does not 
converge on this example, while whitened p-PCR6 usually does. Artificially raising 
the prediction noise solves this problem for ‘standard’ p-PCR6, showing its trend to 


over-concentrate the particle cloud. 


150 


100 






OF True trajectory 
Whitened p-PCR5 ——— 


Particles + 


100 150 200 250 300 


(a) Time step 1. 


150 


100 





OF True trajectory 
Whitened p-PCRS ——— 


Particles + 


100 150 200 250 300 


(c) Time step 20. 


Figure 4.8: The real mobile starts at (200,0) and moves upward at constant 


speed (0,1); bad filters initialization. 
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(b) Time step 10. 


150 


100 


OF True trajectory 
Whitened p-PCR5 ——— 


Particles + 


100 150 200 250 300 


(d) Time step 60. 
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4.4.3.4 Whitened p-PCR6 versus mean 





As seen before, the PCR6-fusion of two probabilistic densities amplifies the areas 
where both densities have a non-negligible value. Otherwise, it usually works like just 
averaging the two densities. In order to measure the impact of the amplification, we 
reprocessed the first example of previous subsection while using the mean, pmean = 
P1t?2 instead of p-PCR6. The result (figure 4.9) is self explanatory: the same 


2 
expansion as with p-PCR6 occurs (figure 4.7), but contraction never happens. 





TEO 


Ley 








DU 
100 150 209 26 300 


(a) Time step 1. (b) Time step 30. (c) Time step 45. 


Figure 4.9: Using mean instead of p-PCR6. Red dots are the positions of the 
particles after fusion. The real mobile starts from (200,0) at time step 0 and 
moves at the constant speed (0,1). 


4.4.3.5 Concluding remarks 


The results presented here have clearly shown that p-PCR6, and especially whitened 
p-PCR6, filters are more robust than the Bayesian filter against bad initialization. 
However, it is clear that Bayesian filters are the best, when the priors are correctly 
defined. The real interest of p-PCR6 is that it does not need a precise prior knowledge 
about the antedating local particle filters. 





4.5 Conclusions 


This paper has investigated a new fusion rule, p-PCR6, for fusing probabilistic den- 
sities. This rule is derived from the PCR6 rule for fusing evidences. It has a simple 
interpretation from a sampling point of view. p-PRC6 has been compared to the 
Bayesian rule on a simple fusion example. Then, it has been shown that p-PCR6 was 
able to maintain multiple hypothesis in the fusion process, by generating multiple 
modes. Thus, more robustness of p-PCR6 were foreseeable in comparison to Bayes’ 
rule. This robustness has been tested successfully on examples of distributed target 
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tracking. It is expected that this new rule will have many applications, in particular 
in case of ill-posed filtering problems. 
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